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1 Introduction 

We will study here the distribution of electromagnetic fields in dielectric periodic media. The 
main difference with the previous topic comes from the word "dielectric" . Obviously even a 2D 
periodic dielectric medium cannot be studied with the Green's functions presented in a previous 
lecture, since the Green's function was for periodic metallic structures. 

In this topic, we will study the EM fields in media where: 

• The material is macroscopic and isotropic. 

• Since the material is constituted of real dielectric, we suppose we work in a small enough 
frequency band such that we can ignore the frequency dispersive behavior of e. 

• The dielectric are lossless, so that e is purely real. 

2 Wave equations 

2.1 Wave equations for H 

Starting from Maxwell's equations and using a permittivity e = 6(f), it is straightforward to 
show that we can write the following equations: 




(la) 



V • H{f) = . 



(lb) 



In this approach, the strategy is therefore: 



1. Find the modes of H{r). 



2. Find those of E{r) by solving Maxwell's equations. 
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2.2 Wave equations for E 



Note that we can rewrite Eq. (la) as an eigenvalue problem: 
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UJ 



@Hif)=[ - ] H{f), 



where 



e = V 



1 



(2) 



(3) 



Upon solving, we get the eigenvectors which correspond to the field patterns of the harmonic 
modes, and the eigenvalues which are proportional to the squared frequencies of these modes. 

Note that O is a linear operator and that it is Hermitian. The demonstration of the last 
property is straightforward: 

Define a scalar product < F, G >= / df F^{f) • G{f) and show that (by 
integration by part for example): < F^QG >=< 6F, (5 >, which is the 
definition of a Hermitian operator. The consequence of this property is of 
course that O has real eigenvalues. 

2.2 Wave equations for E 

Another approach to get the fields would be of course to write the wave equation for E instead 
of for H: 



V X V X E{f) 



e{f)E{f) 



(4) 



However, this system cannot be cast in a simple eigenvalue problem. Although it can still 
be solved for, it is far more complicated to get accurate results since the operator we would have 
to define would not be Hermitian. For this reason, this approach is in general avoided. 



2.3 Bloch states 

In a periodic medium, we know that the fields can be written as: 

Hk (f) = e'^'^Uk (f ) , Uk{f + R) = Uk (f ) . 
Inserting into Eq. (2) yields: 



{ik + V) X 
so that the operator becomes: 



{ik + V) X Uk{r) 



e = {ik + V) 



[e{f) 



(5) 



(6) 



(7) 



Note that because Uk{f + R) = U]^{f)^ the eigenvalue problem can be restricted to a small 
zone in space, which would necessarily imply a discrete spectrum of eigenvalues. Therefore, we 
expect a set of discrete modes for each k. 
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3 Fundamentals of photonic crystals 

We shall briefly explain some terminology here related essentially to solid state physics, but 
which is of prime importance for the study of the structures we are dealing with here. 

3*1 I3irGCt IclttlCG ^some details are given in the "mathematical details" addenda^ 

A photonic crystal is a periodic structure (that we will take to be dielectric here) in ID, 2D or 
3D. 

Any vector f' in space can be written as 

f'^f^R, (8) 

where R is the translational vector in space defined by 

R — aicii + 0:2^2 + <^3tt3 , (9) 

where ai,2,3 ^ {• • • , —2, — 1, 0, 1, 2, 3, . . .} and ai, a2 and as are the lattice vectors. 
From the lattice, we can construct the Wigner-Seitz cell as shown in Fig. 1. 



Figure 1: Wigner-Seitz cell for an arbitrary position of points: the cell is constructed 
by joining the center element to its closest neighbors and drawing perpendicular lines 
from to the center of these segments. The polygon thus created is the smallest repeat- 
able cell of the periodic lattice, and is defined as the Wigner-Seitz cell. 

Note that there exist only one type of lattice for a ID photonic crystal, five distinct types 
for 2D photonic crystals (rectangular, square, hexagonal or triangle, centered rectangular and 
oblique), and fourteen for 3D photonic crystals. 

3*2 r^GCiprOCcll IclttlCG ^some details are given in the "mathematical details" addenda^ 

We will use here the same notation as [Joannopoulos, Meade, and Winn, Photonic Crystals'''] 
and write the reciprocal translational vector as G: 

G = + [32h2 + /Js^a , (10) 

where /3i,2,3 ^ {• • • , —2, —1, 0, 1, 2, 3, . . .} and 61, 62 and 63 are the lattice vectors in the spectral 
domain. 

For the sake of illustration. Tab. 1 gives the definition of vectors a and h for square and 
triangular lattices. 
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3.3 Bloch-Floquet theorem 





Ci\ — dx 
a2 — ay 


U\ — Zj n 1 (Jj X 

62 = 2Ti/ay 


Triangular lattice 


ai — ax 

a2 = (x + V3y) 


hi — 2t{ /a (x — a/3/3 y) 
62 = 27r/a 2^/3 y 



Table 1: Definition of a and h vectors for square and triangular (or hexagonal) lattice. 



3.3 Bloch-Floquet theorem 

From Bloch-Floquet theorem, we know that we can write the electric and magnetic fields as a 
summation over reciprocal vectors G (see the "mathematical details" part of the notes). This 
means that different k do not necessarily correspond to different modes and that therefore there 
is a redundancy in the label k\ we can therefore reduce the study to what is called the first 
Brillouin zone (which is the Wigner-Seitz cell in the reciprocal lattice). Some examples of 
Brillouin zones are given in Fig. 2 and Fig. 3. 

o 



o o 

Figure 2: Direct square lattice and corresponding reciprocal lattice with highlighted 
Brillouin zone. 





Figure 3: Direct triangular (or hexagonal) lattice and corresponding reciprocal lattice 
with highlighted Brillouin zone. 



4 Bragg-like diffraction 

The standard Bragg diffraction is illustrated in Fig. 4. Here, we will derive another diffraction 
condition, equivalent to Bragg, and shall see that the diffraction is entirely governed by the 
reciprocal vector G. 



Figure 4: Schematic representation of Bragg diffraction. Maximal diffraction occurs 
at 2a sin ^ = n\ where A is the wavelength of the electromagnetic wave and n is an 
integer. 

Referring to Fig. 5, we can write the scattering amplitude in terms of the reflection coefficient 
r at position f times a phase factor. 




Figure 5: Diffraction from an elementary volume of a periodic medium: k is the 
wavevector of the incident wave whereas k' is the wavevector of the diffracted wave. 

Upon integrating over the whole volume, we get: 

F{k,k') = j V{f)e'^^-~^'>^dv. (11) 

Since the medium is periodic, we can write: 

r(f + R) = r(f) = f (G')e'^-" , (12) 

G 

such that 

F(k, k')^J2 f t{G)e'^^-^~^>^dv , (13) 

G 
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Section 4. Bragg-like diffraction 



where AA: — k' — k. This amphtude is maximal when G — Ak — 2rmi or, when m = 0, 



(14) 



This is an important relation which, again, is a condition for maximal diffraction. Upon ex- 
panding back in terms of k and k' and rising to the square, we write (noting that = \k'\ = k): 



(15) 



or (taking —k instead of k): 



2k-G = G\ 



(16) 



As an exercise, it is interesting to show that this condition is equivalent to the standard 
Bragg diffraction. 



Upon dividing both terms of Eq. (16) by 4, we eventually write 



(17) 



This last equation has a nice geometrical interpretation shown in Fig. 6 which shows that 
the vectors k that satisfy the maximum diffraction condition are actually those which lie on the 
edge of the Brillouin zone. 



o 



/Gd/2 








O 



Gc/^ ^ C 



O 



O 



o 



Figure 6: Graphical representation of Eq. (17): each vector k (black vector) with its 
tip on a dashed line (not all represented) will satisfy the equation. Graphically: all 
those k have the same projection on the generating vector (5/2 (red vector). 

Therefore: 

The edge of the Brillouin zone plus its center (G = 0) satisfy the maximum 

diffraction condition. 



This condition can also be rewritten in terms of group velocity: for those k which tip lie on 
the edge of the Brillouin zone and = 0, the component of the group velocity normal to the 
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Bragg diffraction planes tends to zero since the electromagnetic wave tends to be completely 
reflected for these k: 

(\ norm 
Vkuj{k)] (fcGBZtip)^O. (18) 

For the symmetry points, the diffracted wave is reflected in the direction of the incident wave so 
that for these points, the total group velocity is zero. This can be directly seen on the dispersion 
curves where, at the symmetry points of the crystal, the tangent to the curve is horizontal 
(except possibly for those points corresponding to a zero frequency). 



5 Mathematical details 

Using all the principles shown before, we can construct the eigenvalue system for H and then 
solve for E. The detailed mathematical manipulations are given in the annex document "Study 
of EM waves in Periodic Structures (mathematical details)". 

Note that to build the system, we need to evaluate the Fourier coefficients of the permittivity 
(or the inverse of the permittivity, k). We shall show how to get these coefficients for the case 
of infinite dielectric rods of circular cross-section organized in a square lattice, embedded in 
a background of e^. We therefore place ourselves in a 2D situation where the G vector will be 
written Gp to denote that it does not depend on z (and similarly R will be noted Rp). 

As a reminder, we write the permittivity as: 

where Q denotes the surface of the elementary cell. 
The idea is to write the permittivity as 

e{p) = et + {ea - e^) ^ S{Rc - \p - Rp\) , (20) 

where again the subscript p in Rp denotes a dependency on x and y only, Rc is the radius of the 
dielectric rods and S denotes the step function. Merging these two equations, we get: 



e{Gp) = ^ f eb + {ea- e^) ^ S{Rc - \p - Rp 



i j efee-^Gp-Pdp + 1 / (e„ - e,) ^ S{R, -\p- Rp\)e-'^^-Pdp . (21) 



Let us call the first integral Ii{Gp) and the second integral l2{Gp)^ and evaluate them separately. 
Evaluation of Ii{p) 
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Section 5. Mathematical details 



The first integral is easily evaluated as: 



€b if Gp = 
elsewhere 



(22) 



Evaluation of l2{p) 

For the second integral, we can make the change of variable p' — p — Rp. Since p spans the 
whole domain and Rp is the translational vector, p^ spans the whole space. We can therefore 
replace the sum of integrals over by a single integral over the whole 2D space. We write then: 



IfGp = 0: 



h{Gp) = II dp'S{Rc -\p- R,\) = (6„ - ^ = friea - et), 



(23) 



where fr = is the fractional volume. 



/2(G'p) = -(€„-€„) 



// dp'S(Rc 

J J — OO 



\p-R^\)e-'^^-p' 



dp'p' / d4>e- 
-/o Jo 



27r j''\p'p'Mp'Gp), 
Jo 



(24) 



where we have used the change of variable (now standard) = p'sin0, = p' coscj)^ 
Gx = GpSin^, Gy = Gp cos 9^ and the well-known identity for the Bessel function. Upon 
using J xJo{ax)dx = x/a Ji(ax), we continue with 



27r 



(Ca - €&) 



G, 



■MGpp') 



Rc 



2Ji{GpRc) 



GpRc 



(25) 



This lead to the final result that: 



~<Gp) = { 



if Gp = 



eafr + eb{l - fr) 

(ca - ^b)fr ^"^G^Rf"^ elsewhere 



(26) 



The reconstruction of the permittivity is then straightforward, and examples of two lattice 
are given in Fig. 7. 
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Square lattice: a = 0.001, e^=10, £^^=1, Rc = 0.0002, (f=0.12566) 




(a) Side view. 

Square lattice: a = 0.001, e^=10, £^^=1, Rc = 0.0002, (f=0.12566) 
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(b) Top view. 

Figure 7: Reconstruction of the permittivity for cylindrical rods in a square lattice 
€a = 10, 65 = 1, Rc/d = 0.2 (yielding a fractional volume fr = 12.5%). 
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Section 6. Dispersion curves 



6 Dispersion curves 

At this point, we have everything to build the eigensystem (34) given in the additional document 
"Study of EM waves in periodic structures (mathematical details)". Solving it gives a set of 
eigenvalues that are directly related to the dispersion curve of the material. An example is given 
in Fig. 8. 




Figure 8: Dispersion diagram for the photonic crystal of Fig. 7. 



